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Abstract. - A recently developed theory of stochastic swimming is used to study the notion of co- 
herence in active systems that couple via hydrodynamic interactions. It is shown that correlations 
between various modes of deformation in stochastic systems play the same role as the relative 
internal phase in deterministic systems. An example is presented where a simple swimmer can 
use these correlations to hunt a non-swimmer by forming a hydrodynamic bound state of tunable 
velocity and equilibrium separation. These results highlight the significance of coherence in the 
collective behavior of nano-scale stochastic swimmers. 



Introduction. 



Swimming strategies for microor- 



^►vganisms and microbots need to take into account the pe- 
r£n culiarities that arise in low Reynolds number hydrody- 
| W-t iamics [TJ-Ej . When utilizing only a small number of de- 
grees of freedom, a careful non-reciprocal prescription of 
*-H cyclic deformations is needed to achieve swimming [1HT5] . 

^"While these ideas have been primarily developed to de- 



[^ 



' scribe the swimming of bacteria |16j , sperms |17j , and 
— * other micro-scale living systems [32] , in recent years they 
£Nj have attracted additional interest with the advent of the 

. * first generation of artificial microswimmer prototypes [19j . 

f- ^ Swimmers of micro- or nanoscale need to face an addi- 
ct t'ional challenge, namely, the overwhelming fluctuations 
^H that would act against their targeted mechanical task. In 
£>. its most basic form, the effect of fluctuations on the mo- 
•'-jtion of swimmers that arc not directionally constrained 
r> or steered is to randomize their orientation via rotational 
^diffusion [201424] . The fluctuations can also interfere with 
the propulsion mechanism and alter the swimming veloc- 
ity [2S], for example via the density fluctuations in the 
case of self-phoretic swimmers [55] or fluctuations in the 
conformational changes in deforming swimmers [27j . 

Interacting swimmers [2 8) are known to have rather 
complex many-body behaviors [29], which can be under- 
stood in terms of instabilities in the context of continuum 
theories (that are constructed based on symmetry consid- 
erations) [30] ■ Another fascinating consequence of long- 
range hydrodynamic interactions between active objects 
with cyclic motions is the significance of internal phase as 



a key dynamical variable [T2"ll3 11152"] , and the possibility of 
synchronization [33]. However, most current theoretical 
studies of the collective behavior of swimmers (continuum 
theories and simulations) ignore the possible effects of co- 
herence, and it is natural to wonder if this is justified. 

One could argue that the overwhelming fluctuations 
that arc present at small scale may wash out any trace 
of coherence among swimmers. To examine the validity 
of this argument, we consider the following question: does 
the notion of relative internal phase apply to stochastic 
swimmers? We use a statistical description to model the 
dynamics of systems that undergo random conformational 
changes while interacting hydrodynamically. Using a spe- 
cific example of a three-sphere system coupled to a two- 
sphere system, we calculate the swimming velocities as 
functions of the statistical transition rates for the confor- 
mational changes. We show that coherence could be intro- 
duced in the system through the correlations between the 
deformations, and that it can be used to create a stable 
bound state between the three-sphere swimmer and the 
two-sphere system (that cannot swim when isolated) with 
a tunable equilibrium distance and velocity. 

Hydrodynamic Model. — Consider the a three- 
sphere system that is located collinearly at a distance 
from a two-sphere system, as shown schematically in Fig. 
[1] The two systems undergo conformational changes by 
opening and closing of the three arms, which could only 
lead to net swimming for the three-sphere system (and 
not the two-sphere system) when isolated, due to scallop 
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theorem [3J. When at a finite distance D, the two sys- 
tems interact hydrodynamically, and their dynamics will 
be coupled to each other. The system on the left (Fig. [T|) 
is made up of three spheres of radius R that are connected 
by arms of lengths L + u^{t) and L + u 2 (t), while the sys- 
tem on the right consists of two similar spheres connected 
by an arm of length L + u R (t). For simplicity, we assume 
that the linker do not interact with the fluid. To Ana- 
lyze the dynamics of the system, we use the linearity of 
the Stokes equation — the equation for hydrodynamics in 
zero Reynolds number — and express the velocity of each 
sphere Vi as a linear combination of the force fj acting on 
a different sphere j: 
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where the details of the hydrodynamic interactions are 
entailed in the coefficients My. Using Oseen's approx- 
imation, we can write simple closed form expressions for 
the coefficients when the spheres are considerably far from 
each other. Denoting the positions of the spheres by #,, 
we have: 

i-nri\xi-Xj\ ' l r Jj 

where 77 is the viscosity of the fluidic medium. Equation 
(fTJ) thus gives us five equation for the ten unknowns Vi 
and /, (i = a, ■■•,e). Maintaining force- free conditions 
on the two systems, namely, f a + ft + f c = and fd + 
f e = 0, provides two additional equations. The final three 
equations are obtained by the kinematic constraints Vb 
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Vb 



•1 • 



and v P 



i\i 



'i R , where the dot 



denotes d/dt. 

Considering the case where the two systems are far from 
each other and the deformations are small compared to the 
average length of the arms L, such that R <C u <C L <C 
D, we can set up a pcrturbative scheme to investigate 
the effect of the hydrodynamic interactions [M] . Solving 
the above linear system of ten equations, we can find all 
the velocities and the forces, from which we can calculate 
the average swimming velocity of the three-bead system 
V L = i (v a + Vb + v c ) and that of the two-bead system 
V R = 2 (vd + v e ). To the leading order in perturbation 
theory, we find 



V L = 



yR 
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We can also extract the average forces acting on the beads. 
To the leading order in deformations, this yields: 

(fa) = \^ (ttf ttf) - 3*^ [(ufu R ) + 2 {u L 2 u R )] , 

(/„> = l^ (tfttf) + 3tt^ [2 (u[u R ) + {u L 2 u R )] , 
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Fig. 1: Schematic view of a three-sphere swimmer interacting 
with a two-sphere system (top) and the three dimensional con- 
figuration space (middle) representing the eight possible dis- 
tinct conformational states of the combined system (bottom). 
We denote by (iot) the state where the three-bead system is in 
state i and the two-bead system is in state a. 



(f d ) = Qvr^ [(uH R ) + (u L 2 u R )] ■ 

Note that (f b ) = -(f a ) - (f c ) and (f e ) = -(f d ). 

The above expressions for the velocities and forces are 
given in terms of the three average quantities (ufu^), 
(u^u R ^, and (u^u 11 ), which correspond to the average 
rates of sweeping enclosed areas in the three perpendicular 
sections of the three dimensional (uf , u 2 , u R ) configura- 
tion space of the system, respectively. For deterministic 
conformational changes of the form 

u\ — dcos(flt — (fi ), 

u 2 = dcos(flt — (p 2 ), 

u R = dcos(D,t- ip R ), 

we can calculate them using time averaging over a period. 
This yields 
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(«M) = ^ 2 Osin(^ - tp%), 

(u[u R ) = -d 2 Qsm(ip R - tp%), 



(5) 

(0) 



Coherent Hydrodynamic Coupling for Stochastic Swimmers 



!£\ 




Fig. 2: Two dimensional projection of the 3D configuration 
space of Fig. [T] showing the different probability current loops. 



(u%u R ) = -d 2 £lsm(<p R - tf%). 



(7) 



The above equations manifestly show that the relative 
importance of these three conformational space area- 
sweeping rates is determined by the relative phases of 
the deformations. We now aim to address the question 
of whether such a concept can exist at small scales where 
the conformational changes are stochastic. 

Stochastic Systems. — To construct a statistical the- 
ory for the deformations of the two systems we assume 
that they have distinct conformational states and the de- 
formations can be modeled as stochastic jumps between 
these states that occur at given rates [37JI33]. The three- 
sphere system can be described with four states and the 
two-sphere system with two states, which make a total of 
eight distinct conformational states in the three dimen- 
sional configuration space, as shown in Fig. [T] More 
specifically the states of the three-sphere swimmer are 
labeled by the index i as follows: i = 1 the two arms 
are closed (uf = 0,u 2 = 0), i = 2 the right arm is 
open (u± = <$*, u% = 0), i = 3 the two arms are open 
(ill = 8,11% = S), and i — 4 the left arm is open 
[u\ = 0,«2 — 8). For the two-sphere system, we only 
have two possibilities: a — 1 the arm is closed (u R = 0) 
and a — 2 it is open (u R = <$*). To describe the instan- 
taneous state of the system we denote the probability of 
finding the left swimmer at state i (i = 1, • • • , 4) and the 
right two-bead system at state a (a = 1, 2) by Pi a . These 
probabilities are normalized as 



X>«» = i. 



The kinetics of the conformational transitions of the two 
coupled systems is given by introducing the corresponding 
transition rates. 



We assume that the conformational changes happen one 
at a time, which means that transitions are only allowed 
between states that are nearest neighbors in the cubic con- 
figuration space shown in Fig. [TJ We denote the tran- 
sition rate for the jump from state i to state j for the 
left swimmer when the two-bead system is in state a by 
fc^(a). Similarly, the transition rate for the two-bead sys- 
tem jumping from state a to state /3 when the three-sphere 
system is in state i is denoted as k R a (i). Note that the 
rates for conformational changes within each system in 
principle depend on the state of the other system. The cu- 
bic configuration space has six current loops corresponding 
to six faces, as shown in Fig. [2] These currents, however, 
are subject to an overall conservation law, which implies 
that only five independent currents exist in the system. 
We define the following currents 

fa = P 11 fc 2 L 1 (l)-P 21 fcf 2 (l), 
fa = fl»l*32(l)--R»l*M(l). 

fa = ftifcii,(i)-.Rii*&(i), 

fa = P n k R (l)-P 12 k R (l), 
fa = Pi 2 fc 2 L 1 (2)-P 22 fcf 2 (2), 

in terms of the probabilities and the rates, and can use 
them to calculate the currents in the loops as follows (see 
Fig. rj: 

h = 77 ( - 3.7 i + J2 + J3 + fa + fa ) , 
6 

h = g (3 ji - 5 fa + fa + fa +fa), 



(3,7*1 + fa -5fa+fa + fa), 



h = g ( - 3.7 1 + .72 + fa - 5 fa + fa ) , 
h = g (3,7*1 + .72 + fa + fa +fa), 



h 



h — I2 — I3 — M — I5 



We can now solve the steady state master equation for the 
system and calculate the probabilities and the currents. 

Using the currents in the loops, we can write down ex- 
pressions for the average rates of sweeping areas in the 
three perpendicular sections of the configuration space, as 
shown in Fig. O The results, which are the statistical 
analogs of eqs. ©, ©, and 0, read 



(u^) = 5 2 (I 6 -I 5 ), 

(uiu R )=S 2 (I 3 -h), 
(u^u R )=5 2 {h-I 2 ). 



(8) 

(9) 
(10) 



Due to the sign convention used in the definition of the 
currents the current running through opposite faces of the 
cube in Fig. [5] have opposite signs. Therefore, the total 
rate of sweeping a certain projected area in the configu- 
ration space is the difference between the currents of the 
corresponding opposite faces in the cube, as eqs. (JSJ), ©, 
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Fig. 3: Phase diagram of possible states. Steady state solutions 
in the form of bound states with constant separation (D = 0) 
are highlighted, and the Stable and Unstable regions are shown. 



and (|T0|) show. The area of each projection is equal to 
S 2 . The above expressions can be used in cqs. (j3|) and (J4]) 
to calculate the average swimming velocities of the two 
systems. 

In the most general case with arbitrary rates the explicit 
form of the resulting probabilities and velocities is cum- 
bersome and therefore not shown here. To illustrate the 
generic features of the solution, we focus on a simplified 
example with the following choices for the transition rates. 
For transitions in the three-sphere system we choose 



kh(a) 



(l + ei)w, i = l,j = 2,a = l, 
(l + e 2 )w, i = l,j = 2,a = 2, 



(11) 



other states 



For the two-sphere system, we choose 

k% a {i) = w, for all states. 



The above choices allow us to only focus on the effect of 
the correlation between the two devices, as having different 
values for t\ and e 2 means that the rate of the three-sphere 
system going from state 1 to state 2, which means opening 
its right arm, depends on whether the two-sphere system 
is in the closed or the open state. If ei = e 2 = 0, detailed 
balance holds and neither of the two components has a net 
motion. With the above choices, the average velocities of 
the two systems can be found as 
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where 



2/0 
2/1 



12(ei + e 2 ) + 5eie2 
56(ei+e 2 ) + 15eie 2 + 192' 

6(e 2 -ei) 

56(e 1 + e 2 ) + 15eie 2 + 192' 



. (12) 
(13) 

(14) 
(15) 



The term proportional to y$ in eq. (|12p is the spontaneous 
swimming velocity of the three-sphere system, while the yo 
contribution in eq. (|13[) is the passive velocity at the loca- 
tion of the two-sphere system caused by the swimming of 
the three-sphere system. The contributions proportional 
to y\ in cqs. (|T2"j) and (TT3")) are active contributions origi- 
nating from a coherence between u\ and u R in the form 
of (uiii R ) ^ 0. Note that the coherence will disappear 
when ei = e 2 . 

An interesting consequence of the coherent coupling is 
that the two systems can form a moving hydrodynamic 
bound state at a fixed separation. Noting that 4fD(t) = 
V R - V L or 

~dT = Ruj {L) 



(2yi - y ) 



.'/o 



(16) 

we can find the conditions at which stable and unstable 
bound states are possible in the effective dynamical equa- 
tion for -D(i), as shown in Fig. [3] The equilibrium dis- 
tance between the two systems in the stable hydrodynamic 
bound states is given as 



D, 



^ 



1/3 



-ei(24 + 5e 2 ) 



12(ei 



£2 



5eie 2 
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(17) 



which can be controlled by changing the transition rates. 
We note that the hydrodynamic-induced formation of 
bound states of a pair of microorganisms has been recently 
observed experimentally |35j . 

Discussion. — Our analysis shows that the concept 
of relative internal phase and coherence between a num- 
ber of systems that undergo stochastic deformations in a 
hydrodynamic medium at low Reynolds number is well 
defined. Stochastic coherence could result from average 
correlations that can be induced between various modes of 
the conformational transitions, and need not exist instan- 
taneously to lead to average correlated behavior. Com- 
parison between eqs. (O, ©, an d ([7]) that are defined for 
deterministic systems and eqs. (JHJ), ©, and (jTUJ) that are 
defined for stochastic systems shows how average relative 
phase between various modes of the conformational tran- 
sitions can be defined and probed in terms of the currents 
in the configuration space of the system. 

In the present study, coherence between the two sub- 
systems is introduced via the rates defined in eq. (fTTl) . 
When ei / e 2 , the rate of opening of the right arm of the 
three-sphere swimmer is chosen to depend on whether the 
two-sphere system is in the closed or open state; it is ex- 
actly this difference that leads to the correlation term pro- 
portional to yi in eqs. (fl"2"|) and (fT3"|) . as can be seen from 
the explicit dependence of y\ ex (e 2 — ei) in eq. (|T5|) . This 
means that while it is possible to have coherence between 
different parts of the system when undergoing stochastic 
fluctuations, this coherence still needs to be imposed via 
the different rates in the kinetic equations. Physically, 
what this means is that the coherence needs to introduced 
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in the system via correlations between different confor- 
mational states of the system and the rates of transitions 
between them. The correlations are reminiscent of the al- 
losteric interactions between proteins [36] . and could in 
principle be engineered for artificial systems using similar 
strategies. These strategies could involve physical interac- 
tions arising from electrostatic forces etc, hydrodynamic 
interactions, and other effects that could modify the tran- 
sition rates by introducing additional mechanical energy 
contributions (costs) in the deformation process and hence 
affecting the transition rates via the exponential (Arrhe- 
nius) dependence on energy change. Alternatively, these 
correlations could be induced via external means such as 
laser pulses that would affect transition rates only in cer- 
tain conformational states. However they are enforced, the 
present study asserts that such correlations could lead to 
a sustainable notion of coherence between stochastically 
fluctuating nano-scale devices in water. 

We have considered the transition rates for the confor- 
mational changes of the two small systems to be time in- 
dependent. If this assumption is not valid for any reason, 
the time dependence in the rates could weaken the degree 
of coherence in the system and ultimately fully eliminate 
it if it is sufficiently strong. This is equivalent to intro- 
ducing time dependence in the phases in the deterministic 
case [described by eqns. (O, ([6]), and 0], which could de- 
stroy the phase coherence. Such a time dependence could 
occur due to temperature fluctuations |37j . which in local 
thermodynamics approximation could affect the transition 
rate through a dependence of the form 



find the following condition: 



u ~ exp 



knT 



(18) 



where / represents a typical force involved in the confor- 
mational change. We can estimate the magnitude of / 
using the typical drag force experienced by a sphere of ra- 
dius R and moving velocity v = Slu, namely / rj 6in]R5 2 u). 
In the present work we have neglected temperature fluc- 
tuations. To examine the validity of this assumption, we 
can use the complementarity relation AE A(l/T) rj — &b, 
which relates the strength of energy and temperature fluc- 
tuations, to estimate the magnitude of temperature fluc- 
tuations as AT 2 = k-gT 2 /C, where C is the heat capacity 
|37j . Using this simplified picture, we can estimate the 
effect of temperature fluctuations on the transition rates 
via 



Alo fS AT 6ttt]R5 2 uj 



k B T 



T 



kvT 



(19) 



In order to have Auj/lo <C 1, which would guarantee that 
the above assumption is valid (i.e. temperature fluctua- 
tions can be ignored) the following condition must hold: 



w < 



fc R T 



(20) 



6irr)RS 2 \ ks 
Putting R = S = 1 nm at room temperature for water, we 



u < (10 8 s" 1 ) x 



(21) 



While there is a debate in the literature about the correct 
choice for C, namely whether it is the heat capacity of the 
entire system [35] or that of the subsystem only [33], eqns. 
(|20p and ([21]) show that in our case neglecting temperature 
fluctuations is justified either way. Note also that our 
analysis has ignored quantum fluctuations, which means 
that the following criterion should also hold: u ^C ^- . 

Finally, we note that we have made the specific choice 
of a three-sphere system and a two-sphere system, because 
its configuration space is 3D and can be easily visualized 
(Figs. [T]and[2]). The same analysis can be easily gener- 
alized to the case of two three-sphere swimmers, in which 
case the configuration space will be 4D and the bookkeep- 
ing of the projected areas in the graph where the proba- 
bility currents are flowing is slightly more delicate. 

In conclusion, we have shown that stochastic swimmers 
could actively couple to each other using hydrodynamic 
interactions. As an example, we demonstrated that the 
coupling could be tailored such that an active swimmer 
could hunt a non-swimming system into a stable moving 
bound state, by "tuning" into the right correlated tran- 
sition rates. Our results show that the notion of internal 
phase and coherence could be important even for fluctu- 
ating systems that are coupled via hydrodynamic interac- 
tions at the nano-scale. 
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